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Abstract 

Background: At the triple border Brazil/Paraguay/Argentina there is easy mobility from one city to another for economic 
and tourism activities. This constant and fast population mobility is mainly to visit Iguazu Falls, in the Iguazu River, on the 
border of the Brazilian state of Parana and the Argentina. As the incidence of tuberculosis is high in this setting, our study 
aimed to establish a first baseline of circulating genotypic lineages of Mycobacterium tuberculosis. 

Methodology/Principal Findings: IHs study included 120 patients from 10 cities in southwestern Parana, Brazil with 
pulmonary symptoms, from July 2009 to July 201 1. Information about sex, age, clinical features and address was collected 
by reviewing the national tuberculosis notification database. Of these, 96 (80%) isolates were identified as M. tuberculosis 
and 22 (22.9%) were drug resistant (20, 20.8% INH mono-resistant and 2, 2.1% multidrug-resistant). All isolates were 
subjected to genotyping by Spoligotyping and MIRU-VNTR typing. The distribution of the isolates analyzed by 
spoligotyping revealed 30 distinct patterns. The four mainly detected clades were Latin American and Mediterranean (LAIVl), 
ill-defined T, Haarlem (H) and S. The MIRU-VNTR showed 85 distinct patterns. Spoligotyping combined to MIRU-VNTR 
allowed 90 distinct patterns. 

Conclusions/Significance: Ou\ study demonstrated that there is significant molecular diversity in circulating M. tuberculosis, 
with predominance of the LAM and T clades in cities of southwestern Parana, Brazil, bordering Argentina and Paraguay. 
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Introduction 

Tuberculosis (TB) is an infectious disease that remains a major 
public health problem in the world. Brazil is among 22 countries 
with major number of cases of TB in the world, with an incidence 
of 46/100,000 populations in 2012 and the distribution of TB 
cases varies gready according to distinct geographic regions [1]. In 
the state of Parana (PR), south of Brazil, 2,350 new TB cases 
(incidence 22.4/100,000) were detected in 2011 [2]. However, 
these rates are higher in certain regions in the state of Parana, such 
as in Paranagua City, a port area, and on the triple border of 
Brazil/Paraguay/ Argentina. The three municipalities concentrat- 
ed near the border, which are Foz do Iguazu in Brazil, Ciudad del 
Este in Paraguay and Puerto Iguazu in Argentina, have an 
incidence of 56.6/100,000, 42.5/100,000, and 23.4/100,000, 
respectively [3]. These high rates are related to the living 
conditions of populations, lack of access to health services in some 



regions and high population mobility. Economic activities and 
tourism to see one of the New Seven Wonders of Nature, Iguazu 
Falls (Cataratas do Iguazu), a waterfall on the Iguazu River on the 
border of the Brazilian state of Parana and Misiones in Argentina, 
have provided strong population growth in Brazilian cities near to 
the border. In this setting, the constant and fast population 
mobility from one city to another can favors the spread of TB. In 
the two major municipalities at the Brazilian border, Foz do 
Igaufu and Cascavel, which have laboratory to perform culture for 
Mycobacterium tuberculosis and attend many others small 
surrounding cities for TB diagnosis had an incidence of TB 
41.8/100,000 and 25.5/100,000 respectively. 

Molecular typing methods have been used for a variety of 
epidemiological investigations and fast tracking geographic distri- 
bution of M. tuberculosis clusters [4]. Identical strains fingerprints 
are called clusters and are usually associated with recent 
transmission, while strains presenting unique fingerprint profile 
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suggest remote transmission [5,6]. A PCR-based method, Myco- 
bacterial Interspersed Repetitive Unit - Variable Number Tandem 
Repeat (MIRU-VNTR) typing allows the high throughput and 
discriminatory analysis of M. tuberculosis clinical isolates [7,8]. 
The 15- and 24-MIRU-VNTR have been indicated for epidemi- 
ological and phylogenetical studies with M. luhnxulosis . However, 
the 12-MIRU-VNTR added of Spoligotyping, another PGR based 
method, has been widely used in the molecular epidemiology of 
TB. Spoligotyping investigates the population structure of M. 
tuberculosis, focusing on the identification of genotypic lineages, 
spoligotype clades and their geographic distribution [9]. Spoligo- 
typing has proven its usefulness for studying phylogeographical 
aspects of prevailing M. tuberculosis genotypic lineages and in 
combination with 12-MIRU-VNTR, it is highly discriminatory 
allowing the study of both the epidemiological and phylogeo- 
graphical aspect of tubercle bacilli, including large-scale studies 
[10,11]. 

This study aimed to establish a first molecular insight of 
circulating genotypic lineages of M. tuberculosis in cities in the 
southwestern Parana State, south of Brazil, which borders the 
countries of Paraguay and Argentina countries, based on 
Spoligot)'ping and 12-MIRU-VNTR typing. 

Results 

Patients, bacterial isolates and drug susceptibility test 

Of 96 M. tuberculosis identified, 64 (66.7'!^)) were isolated from 
male patients, which account for a male/female ratio of about 2:1 
(M/F = 64/32). The age of patients ranged from 18 to 82 years old 
(36.75±13.38) and an age distribution of 17-19 years, 4 (4.1%); 
20-29 years, 25 (26.0%); 30-39 years, 35 (36.7%); 40^9 years, 18 
(18.7%); 50-59 years, 5 (5,2%) and >60 years, 9 (9.3%). Drug 
susceptibihty testing was performed on aU 96 isolates and 
resistance was detected in 22 (22.92%); 20 (20.8%) were INH 
mono-resistant and 2 (2.1'/'o) multidrug-resistant (MDR-TB) 
(Table 1). All resistant M. tuberculosis isolates were from patients 
residing in the municipality of Cascavel, Parana, Brazil. Four of 
these isolates were obtained from prisoners. 

Spoligotyping and MIRU-VNTR typing 

The distribution of 96 M. tuberculosis isolates analyzed by 
spoligotyping revealed 30 distinct patterns (n= 77, 80.2% isolates) 
and 15 (n= 19, 19.8% isolates) that have not been described in 
SITVITWEB database yet. According to SITVITWEB database, 
71/77 (92.2%) isolates were clustered into nine Shared Interna- 
tional Types (SIT), comprising from 2 to 17 isolates each cluster 
and the remaining 6 (7.8%) isolates showed unique SIT. The main 
SITs found in the present study were SIT60 (n=17, 17.7%), 
SIT42 (n=15, 15.6%), SIT53 (n=12, 12.5%), and SIT177 
(n=ll, 11.5%). The 77 (80.2%) isolates, identified in the 
SITVITWEB database, had been classified into clades and 
sublineages, which ranked in four clades, Latin American and 
Mediterranean (LAM) clade and sublineages, ill-defined T clade 
and sublineages, S and Haarlem (H) clades with 50 (52.1%i), 20 
(20.8%), 5 (5.2%) and 2 (2.1%) isolates, respectively (FUe SI). 

Out of 19 (19.8%) isolates none identified in SITVITWEB 
database, 1 1 (1 1.5%) had unique SIT and the remaining 8 (8.3%) 
were included in 4 clusters, comprising 2 isolates for each one 
(Table 1). The clades assignment for these isolates, by SPOT- 
CLUST, revealed tiiat 9 (47.4%) belonged to tiie LAM, 4 (21.1%) 
to T and 3 (15.8%) to H. There were three other isolates 
pertaining to X, S and H37RV (Table 1). 

12 loci-MIRU-VNTR typing conducted in the 96 clinical 
isolates showed a total of 85 distinct MIRU patterns (File SI). 



Seventy-six (79.2%) isolates had unique patterns and the 
remaining 20 (20.8%i) were included in 9 clusters comprising 2 
to 3 isolates for each one (Table 1). Out of 96 isolates, 18 were 
identified in the SITIVITWEB database pertaining to MIRU 
international types (MITs) 19, 25, 33, 39, 42, 140, 161, 201, 218, 
236, 246, 519, 738 and 816 (Tabic 1). 

Allele polymorphism analysis of 12 loci-MIRU-VNTR by 
Hunter-Gaston Discrimination Index (HGDI) [22], revealed that 
MIRU locus 40 was the most discriminatory with 8 alleles, 
followed by MIRU locus 26 with 6 alleles. In the MIRU locus 40 
the presence of single copy and 4 alleles were the most frequent. 
MIRU loci 2, 10, 16, 20, 23, 27 and 31 were moderately 
discriminant. MIRU loci 4, 24 and 39 were the poorly 
discriminant (Table 2). The amplification of MIRU locus 16 of 
isolate 69 was not possible. 

All 96 M. tuberculosis isolates analyzed by spoligotyping and 
12-MIRU-VNTR typing combined could be genotypicaUy differ- 
entiated in 90 patterns. In this analysis, 84 (87.50%) isolates 
showed distinct genotypes and 12 (12.50%) were included into six 
clusters (2 isolates each with 100% similarity) (Table 3). 

Discussion 

This is the first study conducted to establish the circulating 
genotypic lineages of M. tuberculosis isolatc-d from TB patients in 
southwestern Parana, Brazil that borders the countries of 
Paraguay and Argentina. Economic activities and tourism in the 
triple border have provided strong population growth in some 
cities from the Brazil border, mainly in Gascavel and Foz do 
Igua^u. Constant population mobility and the rapid migration for 
commercial purpose, from one city to another one, at the triple 
border, may favor the spread of TB in this setting. 

The age and sex distribution of patients in our study reflects the 
predominance of TB among adult male people (81.4% of all 
patients) as observed in all other regions in the word [1]. The 2:1 
male to female sex-ratio may be explained by the local socio- 
cultural barriers. In this setting, women stay mainly at home and 
men are engaged in external activity, which leads the latter to a 
higher exposition to the bacillus. In the present study, the 
prevalent isolates came from Cascavel, Foz de Iguagu and other 
surrounding small cities. 

The LAM, T and H were the largest clades observed in the 
present study, the three genotypic clades most frequent in Africa, 
Central America, Europe and South America [10]. The present 
results are similar to a previous one, carried out in our laboratory, 
by Noguti et al. [12], which had greater variability in SITs, but 
LAM, H and T were the largest clades in the northwestern 
Parana. In Araraquara, Sao Paulo State, which borders northern 
Parana, Mendes et al. [13] identified 73% of spohgopattems in 
SpoIDB4, and LAM, T and H clades were the most prevalent too. 
More recently, Santos et al. [1 1], in a comparative study between 
1998-2001 and 2002-2006, observed replacement of ill-dcfmcd T 
family for LAM in Araraquara, Sao Paulo. Furthermore, there are 
a number of lineages found in Northern Parana and Sao Paulo 
States studies, including Beijing, which was not detected in the 
present study. 

The most prevalent SITs in southwestern Parana were 60 
(prototype of LAM4), 42 (prototype of LAM9), 53 (prototype of 
Tl) and 177 (prototv^x- of LAM9). These four most prevalent SITs 
observed in our study were found in Argentina too, as reported in 
SITVITWEB and SpolDB4, but no reference to Paraguay was 
made [10]. 

The SITs 334 and 713 (both susceptible and prototype of Tl), 
characterized by one isolate for each one, have not been described 
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Table 3. Clusters by Spoligotyplng/MIRU combination and epidemiological data of 12 Mycobacterium tuberculosis isolates from 
patients in cities in southwestern Parana, Brazil. 



Date of culture/microscopy Epidemiologic 



Isolates 


spoligotyping 


MIRU-VNTR 


of Baciloscopy 


DST 


relationships 


15 


777777777760731 


324326143324 


11/25/2009 


S 


No epimiological relation 


16 


777777777760731 


324326143324 


07/30/2010 


H 


No epimiological relation 


43 


777777607760731 


124326143324 


06/28/2010 


H 


Live in the same neighborhood 
and nearby street 


77 


777777607760731 


124326143324 


11/28/2010 


H 




51 


777777607760731 


124326132324 


07/25/2010 


MDR 


No epimiological relation 


53 


777777607760731 


124326132324 


07/30/2010 


S 


No epimiological relation 


64 


377777607760771 


224236553321 


09/20/2010 


S 


Live in neighboring cities 


66 


377777607760771 


224236553321 


10/20/2010 


S 




85 


377777607760771 


224306153128 


29/12/2010 


S 


No epimiological relation 


86 


377777607760771 


224306153128 


13/01/2011 


S 


No epimiological relation 


107 


776377777760771 


214205163321 


17/05/2010 


S 


No epimiological relation 


120 


776377777760771 


233325153224 


10/05/2011 


S 


No epimiological relation 


DST, Drug 


susceptibility testing; 


H, isoniazid resistant; MDR, isoniazid and rifampicin resistant; S, susceptible to 


soniazid, rifampicin 


and pyrazinamide. 



doi:1 0.1 371 /journal.pone.Ol 071 06.t003 



in Paraguay and Argentina yet. The first one has been described in 
Brazil only once, but the second hasn't been yet (http:/ /www. 
pasteur-guadeIoup.fr:8081/SITVITDemo/). The other spoligo- 
patterns, whicli were described in SITVITWEB database, were 
already identified in Brazil and some of them in Argentina and 
Paraguay. The SITs 52 (prototype of T2), 60 (prototype of LAM4), 
196 (prototype of Tl), 822 (prototype of LAM9) and 866 
(prototype of LAM9) were described for the first time in Parana, 
Brazil (http:/ /www.pasteur-guadeloup.fr:8081/SITVITDemo/). 
Of interest, INH mono-resistant was predominant in isolates 
genotyped as SIT60 (LAM4, 9 isolates), which is frequent in the 
Americas and some countries in Europe (http://www.pasteur- 
guadeloup.fr:8081/SITVITDemo/). The two MDR were geno- 
typed as SIT 60 (LAM4) but had different MIRU patterns. 

The spoligoforest tree indicates that there are conserved 
genotypes as the SIT53 (prototype of Tl), which demonstrated 
the ability to generate 7 other genotypes, by the loss of 1 to 6 
spacers, which predominate T clade. We could observe with 
SIT42 (prototype of IAM9) that the loss of only one spacer, in a 
different section of DR region, yielded 6 different SITs (20, 60, 64, 
177, 822, 866) belonging to the LAM clade (figure 1). 

We observed a total of 9 single locus variation (SLV) events 
within the 12-locus subset. The best discriminatory power was 
obtained considering that a unique SLV in any one of the loci 
analyzed, remains strongly predictive of the absence of a link 
between the isolates. According to Hunter and Gaston [22] the 
most allelic diversity observed in the M. tuberculosis isolates 
studied was in MIRU loci 40 and 26 and moderated polymor- 
phism were found in MIRU loci 2, 10, 16, 20, 23, 27 and 31. 
These allelic diversities also differs from the obtained in 
northwestern Parana by Noguti et al. [12] that observed most 
allelic diversity in MIRU loci 40, 23, 10 and 16 and moderate 
polymorphism in MIRU loci 26, 20 and 31. Sharma et al. [14] 
working with only 6 loci-MIRU (MIRU 4, 10, 16, 26, 39 and 40) 
observed the most discriminatory loci, in order of diversity, 26, 10, 
16 and 40 in isolates from Kanpur, India. Kovalev et al. [15] 
observed most allelic diversity in MIRU loci 26, 31 and 10 in M. 
tuberculosis isolated in Ural region, Russian Federation. On the 
other hand these authors refer to MIRU locus 23 as having the 



lowest discriminatory power, differing from our results, which was 
MIRU 39. 

Despite of high discriminatory power of MIRU typing observed 
in our study, there were isolates clustered by this method. In these 
cases, some clustered isolates (25, 109, 29, 30, 1, 17, 1 1, 13 and 22) 
were discriminated by spoligotyping, demonstrating the need of 
using these two techniques combined to provide a higher 
discriminatory power. In this sense, Spoligotyping combined to 
1 2 loci-MIRU- VNTR allowed 90 different molecular patterns for 
all 96 M. tuberculosis isolates analyzed. In these patterns there 
were 1 2 isolates that weren't possible to be differentiated as unique 
and were included in six patterns with two isolates each one. Four 
these clustered isolates are of particular concern to be discussed. 
The isolates 43 and 77, which were INH mono-resistant, were 




Figure 1. Spoligoforest tree of 96 M. tuberculosis clinical 
isolates from cities in the southwestern Parana, Brazil, showing 
all spoligotyping clades represent changes (loss of spacers) by 
lane. Edges between nodes reflect evolutionary relationships between 
spoligotypes with arrowheads pointing to descendants. The area of 
each node is proportional to the nunnber of isolates. 
doi:1 0.1 371/journal.pone.01 071 06.g001 
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Figure 2. The triple border Brazil/Paraguay/Argentina, with the location of the cities of Foz do Iguacu and Cascavel (Brazil), Ciudad 
del Este (Paraguay) and Puerto iguazu (Argentina). 

doi:1 0.1 371 /journal.pone.01 071 06.g002 



obtained from patients, who lived in the same neighborhood and 
nearby streets and were referred to the central laboratory in 
Cascavel for TB laboratorial diagnosis in the same year. Two 
other isolates (64 and 66) with the same spoligopattern and MIRU, 
were obtained from patients that lived in neighboring cities 
(Cascavel and Campo Bonito) and were referred to the central 
laboratory in Cascavel for TB laboratorial diagnosis. Thus, to 
overcome the observed limitation, additional study should be 
considered by using optimized 15- and 24- loci-MIRU typing to 
differentiate the isolates and provide additional data about any 
isolates relationship. 

In conclusion, our study oflFers the first baseline about the 
genetic diversity of M. tuberculosis isolated from patients with TB 
in cities in southwestern Parana State, Brazil, which borders two 
other countries, Paraguay and Argentina. This kind of study is 
important for tracing relationships among the strains and 



recognition of the mainly clades responsible for spreading the 
disease. One limitation of the present study is the small number of 
isolates from some cities in the samples. However this issue wUl be 
addressed by additional and better designed studies in this setting, 
with a higher timescale and number of isolates conducted in the 
future to understand the chain of TB transmission in this region. 

Materials and Methods 

Patients, bacterial isolates and drug susceptibility test 

The study was approved by the Ethics Committee of the State 
University of Maringa, Parana (protocol No. 018/2011, CAAE 
No. 0375.0.093.000-10). The information about sex, age, clinical 
features and address of patients was collected by reviewing the 
national TB notification database. As the study was conducted 
with M. tuberculosis isolated from sputa that were disposed of, in 
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accordance with internal biosecurity ruli;s, after test results release 
to the patient, and no additional patients' sputa samples and 
information, which could to identify them were requested 
(identification was carried out by a number given to each patient), 
the document participant consents achievement were exempted by 
the ethics committee. The study included 120 patients (66.7% 
male patients, age range from 17 to 82 years old) from 10 cities in 
southwestern Parana, Brazil, which borders Paraguay and 
Argentina countries (Figure 2), with confirmed pulmonary TB 
by microscopy (Ziehl-Neelsen - ZN) and culture (Ogawa medium 
— OK). All patients were attended at two central laboratories in 
Cascavel and Foz do Iguagu cities, southwestern Parana, from July 
2009 to July 2011. Of these, 96 isolates (80%) were identified as 
M. tuberculosis by biochemical methods [16] and drug su.scepti- 
bility tests (DSTs) were carried out for Isoniazid (INH), Rifampicin 
(RMP) and Pyrazinamide (PZA) by Lowenstein-Jensen (L-J) 
Proportion Method [17]. 

Spoligotyping and MIRU-VNTR typing 

M. tuberculosis DNA of 96 M. tuberculosis isolates were 
extracted from cultured cells as described by Cardoso et al. [18] 
and subjected to genotyping by spoligotyping to detect presence or 
absence of 43 spacers by using the standard method [9,19] and by 
the original 12-MIRU-VNTR loci (2, 4, 10, 16, 20, 23, 24, 26, 27, 
31, 39 and 40) as described by Supply et al. [7,20]. MIRU allele 
scoring wa.s determined according to Mazars et al. [21] and 
Supply et al. [8]. The results from each of the 12 loci were 
combined to create 12-digit allelic profiles. The power of 
discriminating 12 loci was determined by Hunter-Gaston 
Discrimination Index (HGDI) [22]. 
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Database analysis 

The spoligo and 12-loci-MIRU patterns were compared to 
international database SITIVITWEB (Institute Pasteur de la 

Guadalupe) [2,3], which is an updated version of the published 
SpolDB4.() database [10] and available at www.pasteur- 
guadeloupe.fr:8081/SITVITDemo/. The spoligopatterns not iden- 
tified in SITVITWEB database were assigned in families by 
SPOTCLUST at http://cgi2.cs.rpi.edu/~bennek/SPOTCLUST. 
html. For additional analysis a phylogenetic tree (spoligoforest) was 
constructed with data from spoligotyping using SpolTools at http:/ / 
www.emi.unsw.edu.au/ spolTools/. 

BioNumerics software (version 4.45; Applied Maths, Srnt- 
Martens-Latem, Belgium) was used for analysis of Spoligotyping 
and MIRU-VNTR patterns. Dendrogram was constructed for 
Spoligotyping and MIRU-VNTR combined. The genetic distance 
was buUt employing the UPGMA algorithm (Unweighted Pair 
Group Method with Arithmetic Mean). 
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